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Abstract 

We present a new general method for proving global decay of energy 
through a suitable spacetime foliation, as well as pointwise decay, starting 
from an integrated local energy decay estimate. The method is quite 
robust, requiring only physical space techniques, and circumvents use of 
multipliers or commutators with weights growing in t. In particular, the 
method applies to a wide class of perturbations of Minkowski space as 
well as to Schwarzschild and Kerr black hole exteriors. 

1 Introduction 

The wave equation 

□g0 = O (1) 

on general Lorentzian background metrics g appears ubiquitously in mathemat- 
ical physics, often in the context of the linearisation of a field theory governed 
by a system of non- linear hyperbolic p.d.e.'s. The most classical example oc- 
curs perhaps in the context of the Euler equations of fluid mechanics, in which 
case g represents the so-called acoustical metric. Another important source for 
(H]) arises from problems in general relativity, where g represents the metric of 
spacetime, and ([T]) can be viewed as a poor-man's linearisation of the Einstein 
equations. For such metrics g, the global causal structure is often much richer 
than that of Minkowski space. In this context, perhaps the most interesting 
examples are so-called "black hole" metrics. Briefly, these are metrics which- 
like Minkowski space-possess a natural asymptotic structure representing future 
null infinity I"^ , but where the past of X"*" has a non-trivial complement in the 
spacetime. 
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A fundamental problem for ([T]) on suitable metrics g is to obtain "decay 
estimates" , that is to say, estimates giving upper bounds for the rates of decay 
for <f) toward the future, in terms of quantities computable initially on a Cauchy 
surface. This classical problem can be studied as an end in itself, but it is 
often an essential step in proving stability results for non-linear equations which 
give rise to <^ as described above. From the latter point of view, there is an 
important trade off between how much decay one proves and how robust is 
the argument to prove it. The history of non- linear wave equations in the 
past 30 years, culminating in part in the monumental stability proofs (see for 
instance [3]) based on robust estimates for equations of the form ([1]), are a 
testament to how important it is to get this balance right. 

The now standard approach to obtaining robust decay estimates proceeds by 
using weighted multiplier estimates (going back to Morawetz [15]) and weighted 
commutators (going back to Klainerman |10j). both originating in one form or 
another from the symmetries of Minkowski space. In Minkowski space, full use 
of this method (applying both weighted multipliers and weighted commutators) 
obtains (a) weighted energy decay, which can be expressed for instance as 

{{t'/r' + 1)0' + {t + r+ l)2((a, + drW 
+ {\t-r\ + ir{{dt-drW + {t' + r')m\')) 

< cf,__,i<i>' + r'iidt^r + {dr^r + m'), (2) 

where denotes angular derivatives, (b) pointwise decay for (/) of the form 

101 < C^{\t -r\ + l)-''Wt + r\ + l)-i (3) 

for t > 0, where E represents a higher order weighted energy at t = 0, and (c) 
pointwise estimates for derivatives of (f>, with additional decay, for instance: 

\dt(P\ < C^/E{\t -r\ + l)-^''-{\t + r\ + (4) 

Obtaining (a) requires only use of inverted time translations {t^ + r^)dt + 2trdr 
as a multiplier. Obtaining the full (b) and (c) requires both use of inverted time 
translations and commutation with the full algebra of commutators associated 
to □. 

Various applications to non-linear problems and other geometric settings 
have required variants of the vector field method where the use of certain vector 
field multipliers and/or commutators is avoided and less information is obtained 
concerning (a), (b), or (c) above. Klainerman and Sideris |11| obtained decay 
results by commuting only with tdt + rdr and angular momentum operators, 
and this, used together with integrated local decay estimates in the style of 
Morawetz |16j . has proven useful in a number of problems concerning obstacles, 
multiple speeds, etc. [Tl[9l[T4]. A common feature of all these approaches, how- 
ever, is the necessity of always applying at least some multipliers or commutators 
with weights growing in t. 
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In the present paper, we give a related but different approach to the problem 
of proving suitable decay estimates for solutions to ([1]) which, as we shall see, 
has in principle a wider range of validity and may be even more useful for non- 
linear problems. In particular, the method here avoids both multipliers and 
commutators with weights in t. Our approach does not start from scratch, 
but begins with certain "more basic" estimates as ingredients, estimates which 
have already been proven (either classically or much more recently) in many 
cases of interest, by various methods. The method by which these ingredient 
estimates are obtained is in fact of no importance for the considerations here. 
As we shall see, our approach applies to the study of ([T]) on Minkowski space 
and sufficiently small perturbations thereof, but also to Schwarzschild and Kerr 
black hole spacetimes, and immediately yields essentially the decay consistent 
with non-linear stability. The common feature which these regimes share are 

(i) An integrated local energy decay estimate 

(ii) Good asymptotics towards null infinity allowing for a hierarchy of r- weighted 
energy identities 

(iii) Uniform boundcdness of energy 

Estimates of the form (P, already referred to above, were first proven by Mo- 
rawetz |16| for perturbations of Minkowski space, while for Schwarzschild, es- 
timates of the form Q are due to several authors (see Section 4.4.1 of [5]). 
For slowly rotating Kerr spacetimes {M,gM,a) with parameters \a\ « M, such 
estimates have been shown by [H [17] (see also [2]), whereas, for the general 
sub-extremal Kerr case \a\ < M , this will appear in a forthcoming paper. In 
the black hole case, the "integrated decay" estimate (P, to be useful here, must 
be understood in the sense of an estimate which does not degenerate at the 
horizon T^"*", an estimate which is indeed possible in view of a multiplier con- 
struction intimately connected to the celebrated red-shift effect (see [HIl]). In 
the Schwarzschild and Kerr case, the integrated decay estimate necessarily de- 
generates, however, near the analogue of the photon sphere, and this loss is 
related to the trapping phenomenon; a non-degenerate estimate can neverthe- 
less be obtained after losing differentiability. Similar phenomena occur in the 
simpler setting of obstacle problems in Minkowski space, to which this method 
also applies. 

We note that the hierarchy of r-weighted energy identities (ii) applied here 
corresponds to the use of a multipHer of the form rP{dt + dr) for < p < 2. 

In principle, our approach can retrieve the "full decay statements" of the 
traditional vector field method as would follow from the combined use of [T3] 
and [To] , except that the energy decay statements will in general lose derivatives 
if the integrated local energy decay statement does so. In this short note, we 
shall here only give the first step of the approach. This begins by obtaining 
T^^ decay of the energy flux through a suitable foliation Er (whose leaves ap- 
proach null infinity), where r is the natural parameter of the foliation. This is 
precisely the analogue of (See the discussion of Section[7J) This then leads 
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immediately to pointwise results (in Section [5]) which are the analogue of 

\r^'^^\<Ct-\ \r(t>\ < Ct-^/\ 

Again, no use is made of commutation with vector fields with weights in t. 
We outline in Section [8] how to extend to the method to achieve the analogue 
of ®, &. 

We consider the case of Minkowski space first in Section [51 followed by 
Schwarzschild in Section [H to illustrate the structure of the argument in each 
of these settings. As shall be apparent, the nature of the argument is such 
that elements (p])- (pll)) . suitably interpreted, are sufficient and a very general 
theorem can in fact be formulated. This general theorem includes in particular 
the Kerr case. See Section [HI We leave the precise formulation of the most 
general assumptions for a follow-up paper. 

2 Notation 

Let {Ai ,g) be a Lorentzian manifold and (j) a solution to ([T]) on Ai . We recall 
the energy-momentum tensor 

Given a vector field X , we define the currents 
where tt^^, = ^Cxgtj.iy- Recall that 

If X is Killing, then iir"^[(/)] = 0, in which case is a conserved current. 

Integrands without an explicit measure of integration are to be understood 
with respect to the induced volume form. For the case of null hypersurfaces S, 
integrands of the form J^nf^ are to be understood with respect to a choice 
of volume form on E and corresponding normalisation for such that the 
integral represents the correct term in the energy identity. Alternatively, one 
can interpret this simply as the integral of the 3- form f^*Jx, where * denotes 
the Hodge star operator associated to the spacetime metric. 

3 The Minkowski case 

Let us start by considering the wave equation in Minkowski space 

□ ^=0. 

For convenience, we may assume for this discussion that (j) is smooth and com- 
pactly supported on a Cauchy hypersurface. This will ensure that all integrals 
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are a priori defined. By a density argument, our results will only require the 
finiteness of the quantity on the right hand side of the final estimate ([8]) and 
the relation ^ at spacelike infinity. 

Let {t,r,6,(p) denote a choice of standard spherical polar coordinates on 
Minkowski space. We shall let 'f denote the induced covariant derivative on the 
spheres of constant r, and ^ the induced Laplacian. We will also define the null 
coordinates u = t-r, v = t + r. Finally, let us introduce the notation T to denote 
the coordinate vector field dt with respect to {t, r) coordinates. 

For our argument it will be useful to consider the following foliation: Fix 
R > and consider the hypersurfaces Et defined as the union of the spacelike 
{t = T}n{r<R} and the null {u= {T-R)}n{v> (t + R)}. 

The following estimate can be shown by the original method of Morawetz [16] : 

f°° f Cfl f Jj[0]ng^. (ILED-Mink) 

Jt Jr<R JS-r 

The initials of the label stand for "integrated local energy decay" . Since = T" 
in r < i?, we could replace the first term on the right hand side more explicitly 
with the quantity (dfcj))^ + {dr4))^ + ll^0Pj but we prefer this more geometric 
formulation to compare with other situations. Our goal is to use pLED-Mink[) 
to derive decay for the energy flux through I],-, and finally, pointwise decay 
estimates. 

We rewrite the wave equation in the form 

-dudyip + ^ip = 0, i' ■■= rcj), 

multiply the equation by r^dyip and integrate by parts in the region ZJJJ bounded 
by the two null hypersurfaces ui = (n - R), U2 = {t2 - R) and the timelike 
hypersurface r = R: 



We obtain 



J rP{dyi;f sine d0d(l)dv 

U=T2-B,V>T2+R 

+ f rP-^p{dyijf + {2- p)\fTpf) sine ded(l>dudv+ J rP\fxljf sine de d(l>du 



J rP{dy^l^f sine de del) dv 



u=ri-R,v>Ti+R 

T2 

+ 

Tl 



> 2 

f rP {\f^\^ - (dy^f) sine ded(j)dT\r=R . 



(p-WE-Mink) 



The initials "WE" denote "weighted energy". We have written exphcitly the 
measure of integration to emphasize that the expected part of the volume 
form is now in fact included in the quadratic terms in (remember; -0 = r0). 
Thus, the factors r^, r^"^ are to be viewed as weights. We see that the left hand 
side of the above identity is positive definite for p<2. 

We note that identity ( |p-WE-Mink[ ) can be reinterpreted in terms of the 
energy identity in satisfied by the currents [</>], if^[(/>] corresponding to 
V = rPdy (where dy is understood with respect to coordinates), suitably 

modified, however, by appropriate O'th order terms. 

Let us first apply ( p-WE-Mink[) with p = 2. Then, observing that the last 
term on the right hand side of (p-WE-Mink) can be controlled (after a bit of 
averaging in R) by the left hand side of (jlLED-Minkp , we obtain 

J r^idytpf sine d9d(j)dv+ J r{dy^pf sine d9 d(t)dudv 
f {dytpf sine ded(j)dv + C f 



u=T2-R,v>r2+R 
< 



(5) 



U=Ti -R.V>Ti+R 



This implies that we can find a dyadic sequence of r„ ^ oo with the property 
that 



/ 



r{dvip)^ sin 9 de d(f> dv 



U=T„-R,V>Tn+R 



<Ct- 



J r^{dyil)f sine d9d(j)dv + Jl[(j)]i 



U=Ti -R,V>Ti+R 
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We now apply ( |p-WE-Mink ) with p = 1 to the region to obtain 

J r{dyi}f sine ded(j)dv+ J {{dyi^f + sin 9 d9 d(t)dudv 



U=Tn-R,V>rn + R 



Dl 



f r{d^,i}}fsineded(t)dv + C f jj[</)]rig 

J r^idyijf sine ded(j)dv + J jj[0]n^ 

l.V>Ti + 



U=rn-l~B,,V>Tn +R 



<Ct, 



U=Ti -i?,,lJ>Ti +i? 

Adding a multiple of the estimate pLED-Mink[) . we obtain 

JTn-l Jr<R 

J r{d^'ipf sine ded(l)dv+ j {{dyijjf 



') sine de d(j)dudv 



U=Tji -R^V>Tji +R 



<Ct- 



J r\dyi;f sine ded(l)dv + ^l[0]"s. 



U=Ti -R,V>Ti+R 

a 



(6) 



We now observe that 



/ {d^{r(l))y dv = / [r^{d^(f))^ + 2rd^(P<j>+(t)^]dv= / r^{d^(l)fdv-r(l)'^\^=^„ 
Substituting this into ^ we obtain 



f r'^{d^i}}f sine ded(t)dv+ f 

(7) 



U=Tl -R^V>Tl+R 



Finally, in view also of the energy bound 
which holds for all t > t' , the estimate ([7]) easily implies that 



(EB-Mink) 



for any t > ti . 



<Cr- 



U=Ti ~R,V>Ti +R 



(8) 
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Remark 1. Note that (0) essentially corresponds to (0) of the discussion of 
the vector field method. In view of the absense of weights in t in the multi- 
plier used for the proof it is already immediate from the above that (0) follows 
similarly on sufficiently small non- stationary perturbations of Minkowski space. 
This estimate is in itself very useful for applications to global existence problems 
for quasilinear equations, and previously was not available without commutation 
with inverted time translations with their quadratic weights in t. See the discus- 
sion in Section^ 

4 The Schwarzschild case 

We now consider the wave equation ([T]) on a Schwarzschild background ( , g) 
with parameter M. See Section 2 of [8]. Let V denote the closure of (a connected 
component of) the domain of outer communications, let S,- be a translation- 
invariant family of hypersurfaces connecting the event horizon 'W' and null 
infinity I"*" as below: 




For definiteness, one can choose the family Sr explicitly as follows: Given a 
constant R > 3Af, first define a Regge- Wheeler coordinate r* = r + 2A/log(r - 
2M) - i?, the null coordinates u = t-r*,v = t + r* and the coordinate t* = t + 
2Af log(r-2Af). Then can be defined to coincide with t* = t for r < R, and to 
coincide with u = t for r > R. Let us define -D^J analogously as before, i.e., with 
the present normalisation of the null coordinates: D^J = {r > R} n {n < u < T2}. 

It will be useful to fix a translation-invariant timelike vector field N on 
J'^(I]t-i) such that N = dt for r > R, say. (The coordinate vector field dt* in 
(r, <*) coordinates will do for instance.) Let us note that (jEB-MinkI) is replaced 

by 

Ja[<P]n^. <cf^^ J^i^n^^,. (EB-Schw) 

This non-degenerate energy boundedness statement was originally proven in [5] 
but can in fact be shown independently of decay results, as in U and Section 3 

of 0. 

The analogue of (jILED-Minkp is 

f°° f {xJa[^]n%^ + (drcjyf + 4>^)<C f J^[0]ng^, (ILEDl-Schw) 
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where x is a weight which vanishes quadratically at r = 3M. This type of 
estimate, which does not degenerate on the horizon, was originally proven in [5]. 
For background, see Section 4.4.1 of [5]. Commuting with T, we may obtain for 
instance 

f J^[cl>]nt<C f J^[mn%^+C f J^[</>]ng^. (ILED2-Schw) 

Jt Jr<R J2lr Jllr 

Wc will outline more explicitly further on in this section how the above 
are used in the Schwarzschild case. Let us immediately remark, however, 
that (jILEDl-SchwP would be sufficient for the use of (jILED-Minkp to obtain 
the analogue of (l5|) , whereas one would have to use (jILED2-SchwP to obtain ([6]) , 
necessitating the appearance of a higher order quantity on the right hand side. 
This will mean that, in the Schwarzschild case, the final estimate loses differen- 
tiability, as expectedQ In this latter use of (jILED2-Schwp . it is essential that 
one has the non-degenerate quantity J^l^ffl on the horizon. 

Before adapting the argument of Section [3] to Schwarzschild, it remains to 
derive the analogue of ( p-WE-Mink ) . 

As before, we write the wave equation ([T]) as 



-dud^il) + (1 



2M 



2M(l-^) 



^ = 0, 



if) := r(j). 



We may rewrite this in the form 

1 



1- 



2M 



2M 



-Tp = 0. 



Observe that this implies that 
1 



-du 



1 



2M 



2MrP 



1„ /2MrP\ ,2 
'2^A-^V 



■ rP/^i/jd^i/j = 0, 



which generates the following additional terms in the analogue of ( |p-WE-Mink[ ) : 
a boundary term at null infinity of the right sign, proportional to 



2M 



I 



' sin dO d6 du 



and the bulk terms 



M J dy [r''-^]i;^ sin 9 d9d<f)dudv = {3 -p)M J rP'^^p"^ sin 9 



DZ 



Idvipl"^ sin 6 d6 d(l)dudv ■ 



d9 d(f> du dv 



Idyip]'^ sin9d9 dtp dudv. 



^The precise loss of differentiability in (IILED2-Schwl l is wasteful. Using a refined estimate 
(see 1131 ) would lead to an improved result here in the sense of loss of differentiability, but 
this would take us outside the realm of physical space methods. 
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Both expressions have the right sign as long as < p < 3 and r is sufficiently 



large. Note also the new (1 - 2M/r) weights in the analogue of (p-WE-Mink), 
which of course will play no role for large r. 

Suppressing the spherical integration from the notation, wc thus obtain, for 
R sufficiently large, the estimate 

J rP{d^,il^fdv+ J rP^^{p{d^,i;f + {2-p)\fi)f)dudv+ J 

U=T2,V>T2+R 



<c 



f rP{d^4'fdv + J rP{\t^\^ + {d^il^f)dT\r-.R 

) 

(p-WE-Schw) 



Given (jlLEDl-Schwp . pLED2-Schwp . ( |p-WE-Schw] ), (|EB-SchwD . one repeats 
the Minkowski argument until before 

At this point, one adds the estimate (jILED2-Schw[) to obtain 

T" f +<^')+ f r{d,4'fdv+ f {{d,,i^)^ + \M'')dudv 

JTn-i Jr<R J J 



U=Tn~R,V>Tn+R 



< CtZ^ 



U=Ti -R,V>Ti-\-R 



(9) 



One obtains as before 



U=Ti -R,v>Ti-\-R 

'''n-l '''n-l 

Identifying a good dyadic sequence, applying the boundedness statement (jEB-Schwp . 
and using the above inequality again with replaced by T(f>, we obtain finally 

fj^[<l^]nt<Cr-'j:f rV„^[T>K^. (10) 

5 Pointwise estimates 



Starting from the energy decay bounds (|8]), (|10p proven above, one can obtain 
pointwise decay estimates after additional commutations and applications of 
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weighted Sobolev inequalities. To obtain decay near null infinity, we shall require 
commutator vector fields with weights in r, but, again in accordance with the 
philosophy of our method, no weights in t will be used. 

In the Schwarzschild and Minkowski space, one approach for this is to use 
only commutation with the so-called angular momentum operators fti, defined 
as a standard basis for the Lie algebra corresponding to rotations. See [S]. 

An alternative argument (see Section 5.3.8 of [8]) uses commutation with 
cut-off versions of fli, supported only near infinity, together with additional 
commutations with T and the red-shift vector field N (or equivalently Y) and 
an appeal to elliptic estimates. This argument is more robust, and moreover, can 
also be used to obtain non-degenerate pointwise estimates of higher derivatives 
up to and including the horizon. 

From this latter approach, we obtain finally: 

\r^'^(l>\<C\/ET-\ \r(j)\<C\/ET-^l^, (11) 

where 




We note that the second two terms of the integrand arise from the loss in 
(jILED2-SchwP and are thus unnecessary in the case of Minkowski space, where, 
in addition, of course, we take N = T. 

6 A general formulation and the Kerr case 

It should be clear from the structure of the above argument that the properties 
(i), (ii), (iii), as captured by suitable estimates in the form of (jlLEDl-SchwP - 
(|ILED2-Schwp . dp-WE-Schwl , (jEB-Schw^ . are in fact suflicient for obtaining 
pointwisc-in-time energy decay of the form (llOp . In particular, we note that 
the stationarity of the metric is not used directly, as long as one has the above 
estimates (i), (ii), (iii), and (in the case where there is trapping), as long as 
one can absorb the errors arising from commuting with T. We shall give a 
general formulation in a follow-up paper. The results apply to a wide class of 
background metrics and can include obstacles with trapped rays, black holes, 
etc. 

Of particular importance are the Kerr spacetimes {M,ga,Ai)- See Section 5.1 
of [8] for a discussion. In the case \a\ « M, the estimate associated to (ii) follows 
by the results of [H [17] , whereas the estimate associated to (iii) follows from [7] . 
For the general case \a\ < M, the estimates associated to (i) and (ii) will follow 
by forthcoming results. On the other hand, a computation reveals that the 
analogue of (i) follows as in the Schwarzschild case. Thus, one obtains precisely 
()10p for Kerr, where T is the stationary Killing field, is a globally timelike 
translation-invariant vector field, and fli are angular momentum operators on 
an ambient background Schwarzschild metric. 
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In the general formulation, to obtain pointwise estimates as in (jlip . the 
assumptions necessary for (|10p have to be supplemented with assumptions re- 
garding the existence of sufRciently many good commutators. Again, we leave 
the most general formulation to a follow-up paper. Wc remark here, however, 
that the necessary properties can be immediately explicitly verified for Kerr, for 
the entire range \a\ < M, in view of the timelike nature of the span of dt and 

outside the horizon, and the positivity of the surface gravity ensuring the 
validity of the commutation Theorem 7.2 of [5]. Thus, as in Schwarzschild, one 
obtains the pointwise (fTT|) for Kerr, where the vector fields are interpreted as 
in the previous paragraph. 



7 Discussion 

As mentioned before, the energy decay estimate ([5]) corresponds exactly to what 
one can obtain in Minkowski space by applying the multiplier associated to 
the generator of inverted time translations 

Z = udu + v^dv = {t^ + r'^)dt + 2trdr. 

This argument for obtaining decay estimates in Minkowski space goes back to 
Morawetz [T3| . Analogues of the current can be constructed in the Schwarz- 
schild case [SI E] and slowly rotating Kerr case |a| « M [H [2], and were used in 
conjunction with an integrated decay type inequality, to obtain decay of energy 
in St-. Already, however, in the Kerr case, there is a loss in the power of decay 
which one can obtain by this method, which is in principle proportional to the 
parameter a. See Section 5.3.6 of [8]. The difficulty with these constructions 
arises from error terms in the region r < R which inherit large weights in t. These 
problems could be even more profound in the setting of non-linear applications. 
This issue was in fact our motivation for the present work. 

It is amusing to compare the role of null infinity in the present argument 
with the role of the cosmological horizon in the Schwarzschild-de Sitter case 
studied in [6]. There, the analogue of the iteration process applied here could 
be continued indefinitely to yield decay bounds faster than any polynomial rate. 
The reason is that in the analogue of ( |p-WE-Schw[ ) , the weights of the boundary 
terms coincide with the weights of the bulk terms. This is precisely related to the 
red-shift. Indeed, this is reflected already in the Schwarzschild case at the event 
horizon by the fact that J^[(/>] appears without degeneration at r = 2M on both 
sides of inequality (jILEDl-Schw|) . A similar relation holds at the cosmological 
horizon of Schwarzschild-dc Sitter. In contrast, in the present asymptotically 
flat case, the weights in ( |p-WE-Schw[ ) arc polynomial in r and drop by 1 power 
when comparing bulk and boundary terms. This is why it is important to keep 
the whole p- hierarchy of estimates ( |p-WE-Schw ) in mind for the method to be 
useful. 

Wc mention finally a forthcoming Fourier-based method of Tataru and col- 
laborators for obtaining pointwise estimates for stationary metrics starting again 
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from (i) and good asymptotic behaviour of the metric^. 



8 Refinements: the full decay of the vector field 
method 

As explained in the introduction, the "full vector field method" yields faster 
pointwise decay of the form ([3]) . One way to prove such a result starting from our 
(fTU|) would be to apply one commutation with the scaling vector field tdt + rdr ■ 
(See work of Luk [12] for the Schwarzschild case.) This would be contrary, how- 
ever, to the philosophy of our current approach, which is to avoid any multipliers 
and commutators with weights in t. It turns out that our approach can in fact 
be extended to yield more decay without such commutations, by instead com- 
muting with dy (again understood in a null {u,v) coordinate system), and Qi 
which allows applying the hierarchy of estimates ( |p-WE-Mink I or (p-WE-Schw) 



for larger weights in r, i.e. for larger values of the parameter p. This leads to 
additional decay for a higher-order energy fiux through E,-, which can then be 
transformed into the usual decay for (j) using also ideas of [12]. This will be 
discussed in a forthcoming follow-up paper. 
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